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ABSTRACT 

We report the identification of a recurrent ultraluminous X-ray source (ULX), a highly absorbed X-ray source 
(possibly a background AGN), and a young supernova remnant near the center of the starburst galaxy M82. 
From a series of Chandra observations taken from 1999 to 2005, we found that the transient ULX first appeared 
in 1999 October. The source turned off in 2000 January, but later reappeared and has been active since then. 
The X-ray luminosity of this source varies from below the detection level (~ 2.5 x 10 38 ergs s" 1 ) to its active 
state in between ~7 x 10 39 ergs s _I and 1 .3 x 10 40 ergs s _I (in the 0.5-10 keV energy band) and shows unusual 
spectral changes. The X-ray spectra of some Chandra observations are best fitted with an absorbed power- 
law model with photon index ranging from 1.3 to 1.7. These spectra are similar to those of Galactic black 
hole binary candidates seen in the low/hard state except that a very hard spectrum was seen in one of the 
observations. By comparing with near infrared images taken with the Hubble Space Telescope, the ULX is 
found to be located within a young star cluster. Radio imaging indicates that it is associated with a H II region. 
We suggest that the ULX is likely to be a > 100M© intermediate-mass black hole in the low/hard state. In 
addition to the transient ULX, we also found a highly absorbed hard X-ray source which is likely to be an 
AGN and an ultraluminous X-ray emitting young supernova remnant which may be related to a 100-year old 
gamma-ray burst event, within 2 arcsec of the transient ULX. 

Subject headings: black hole physics — galaxies: individual (M82) — supernova remnants — X-rays: binaries 
— X-rays: galaxies 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are defined as off- 
nuclear X-ray sources with isotropic luminosities much 
higher than the Eddington limit for a solar mass black hole 
(Lx ~ 1.3 x 10 38 ergs s" 1 ). Typical X-ray luminosities of 
ULXs are in between 10 39 ergss _1 and 10 41 ergss _I . The 
physical nature of ULXs has been an enigma because of 
their high energy output. Many ULXs show strong variabil- 
ity suggesting that they are accreting compact objects. As- 
suming the emission is isotropic, then some of the ULXs 
may harbor intermediate-mass black hole (IMBH; Colbert & 
Mushotzky 1999; Makishima et al. 2000) with masses of 
100- 10000M Q . Alternatively, ULXs may simply be stellar- 
mass black holes. It has been suggested that ULXs are stellar- 
mass black holes with radiation pressure-dominated (Begel- 
man 2002) or slim (Ebisawa et al. 2003) accretion disks 
that cause super-Eddington luminosities. Furthermore, ULXs 
may be stellar-mass black holes with anisotropic X-ray emis- 
sion (King et al. 2001), or micro-blazars which happened 
to be observed along the direction of their relativistically 
beamed jet (Kording et al. 2002). In addition, some ULXs 
may be young X-ray luminous supernova remnants in a high- 
density medium, or hypernova remnants. Finally, some ULXs 
have been identified with background AGNs through optical 
follow-up spectroscopy (e.g., Foschini et al. 2002; Masetti et 
al. 2003). Each of these models has difficulties to fully ex- 
plain the observations, but yet has some supporting pieces of 
evidence. Currently, we do not have a complete picture about 
the physical nature of ULXs, primarily because we do not 
have dynamical mass measurements of the compact objects 
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that power ULXs. 

In our current stellar formation and evolution theory, we 
are only able to constrain two classes of black holes, the su- 
permassive black holes with masses exceeding 10 6 M Q at the 
center of galaxies and stellar-mass black holes with masses 
lower than 20 M Q . If some ULXs host IMBHs, they might 
provide a clue to fill in the missing link between stellar-mass 
black holes and supermassive black holes. 

In this paper, we report on a recurrent transient ULX in the 
starburst galaxy M82 by using archival Chandra and Hub- 
ble Space Telescope (HST) data. The source is located near 
the galactic dynamical center and is one of the most lumi- 
nous X-ray sources within the supper-bubble region of M82 
(Matsushita et al. 2005; Matsumoto et al. 2000), close to sev- 
eral super-star clusters. This source is probably the second 
most luminous source in M82. The most luminous source, 
M82 X— 1, is a prime candidate of IMBH and is about 5 arc- 
sec from the transient (see Kaaret et al. 2006 and references 
therein). In addition, we also study the physical nature of the 
two ULXs very close to the ultraluminous transient by using 
multi-wavelength data. 

In § 2 we describe the Chandra and HST observations. We 
present the data analysis and results in § 3. A discussion of 
the three ULXs is given in § 4. 

2. OBSERVATIONS 
2.1. Chandra 

M82 (NGC 3034) is a nearby starburst galaxy. We adopt a 
distance of 3.6 Mpc to M82 based on the Cepheid distance of 
3.63 ±0.34 Mpc to its close neighbor galaxy M81 (Freedman 
etal. 1994). M82 was observed twelve times between the year 
of 1999 and 2005 with Chandra. The details of the observa- 
tions are given in Table 1 . Among these twelve observations, 
observations 3 and 5 were using the High Resolution Camera 
(HRC-I). The rest were using the Advanced CCD Imaging 
Spectrometer array (ACIS-I or ACIS-S). We used CIAO v3.3, 
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TABLE 1 
Chandra OBSERVATION LOG 



Index 


Dale 


Obsld 


Exposure 


Instrument 


Remark 


1 


1999-09-20 


361 


33.7 ks 


ACIS-I 




2 


1999-09-20 


1302 


15.7 ks 


ACIS-I 




3 


1999-10-28 


1411-1 


36.3 ks 


HRC-I 




4 


1999- 12-30 


378 


4.2 ks 


ACIS-I 


off-axis 


5 


2000-01-20 


1411-2 


17.8 ks 


HRC-I 




6 


2000-03- 11 


379 


9.1 ks 


ACIS-I 


off-axis 


7 


2000-05-07 


380- 1 


3.9 ks 


ACIS-I 


ofr-axis 


8 


2000-06- 12 


380-2 


1.2 ks 


ACIS-I 


off-axis 


9 


2002-06- 18 


2933 


18.3 ks 


ACIS-S 




10 


2005-02-04 


6097 


58.2 ks 


ACIS-S 


off-axis, 1/8 subarray 


11 


2005-08- 17 


5644 


75.1 ks 


ACIS-S 


1/8 subarray 


12 


2005-08- 18 


6361 


19.2 ks 


ACIS-S 


1/8 subarray 



NOTE. — ObsIDs 1411 and 380 have two observations merged in one 
event list. We used a time filter to separate the two observations. 

HEAsoft v6.2, and XSPEC vll.3 packages to perform data 
reduction and analysis. 

For ObsIDs 1411 and 380, there are two separate observa- 
tions within the same event list. Therefore, we used a time- 
filter to split the observations and analyzed the data separately 
(observations 3, 5, 7, and 8). Five of the observations (4, 6, 7, 
8, and 9) are off-axis to reduce pile-up of M82 X-l due to its 
high luminosity. In particular, during observations 10-12, the 
detector employed a 1/8 subarray mode with a frame time of 
0.441 s to reduce pile-up. 

2.2. HST 

Since M82 has high extinction near its core where our tar- 
gets are located, we used near IR image to study the envi- 
ronment around the X-ray sources. We obtained HST Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) 
data from the Multimission Archive at STScI (MAST). M82 
was observed on 1998 April 11 by using the NIC2 camera 
with the F160W filter (//-band). The NIC2 camera has a res- 
olution of 0".075 pixel" 1 and a field-of-view of 19.2 x 19.2 
arcsec 2 . To cover the central 1.5 arcmin along the semi- 
major axis, we obtained 1 1 pipeline processed NICMOS im- 
ages from the MAST. The individual images were combined 
to form a mosaic that were used for analysis. To correct the 
absolute astrometry of the NICMOS image, we aligned the 
NICMOS image with a wide-field HST ACS mosaic image 
(Mutchler et al. 2005). The ACS mosaic has a dimension 
of 10'. 24 x 10'. 24 and can be obtained as High-Level Science 
Products via the MAST 4 . The observations and data reduction 
are described in Mutchler et al. (2005). We used the F814W 
(/ band) mosaic for our analysis and corrected the absolute as- 
trometry by using the 2MASS catalog. We identified 21 iso- 
lated stars in the field and matched with the 2MASS catalog. 
Using IRAF tool ccmap, we corrected the absolute astrome- 
try of the ACS mosaic with a registration error of 0.2 arcsec. 
We then registered the NICMOS mosaic with the astrometric 
corrected wide-field ACS image. With 9 isolated stars in both 
field-of-views of NICMOS and ACS mosaics, we corrected 
the astrometry of the NICMOS image with a registration er- 
ror of 0.036 arcsec. 

3. DATA ANALYSIS AND RESULTS 

3.1. X-ray Imaging 

The brightest source in the field of M82 is M82 X-l 
(CXOU J095550.2+694047; see Fig. 1). We note that 
the X-ray coordinates are based on an astrometric corrected 

4 http://archive. stsci.edu/prepds/m82 



image (see § 3.5). About 5 arcsec southeast of M82 X- 
1, there is a complex of three bright X-ray sources (Mat- 
sumoto et al. 2001). By examining the Chandra images, the 
brightest one, CXOU J09555 1.0+694045 (J095551.0 here- 
after; note that in Matsumoto et al. 2001, the source is called 
J09555 1.1 +694045 due to different astrometry), was clearly 
below the detection limit in two observations indicating that 
it is a highly variable source. Note that previous studies 
have mentioned the transient behavior of this source (Mat- 
sumoto et al. 2001; Kaaret et al. 2006). In addition to 
J095551.0, there are several additional transients as shown 
in Figure 1; the discussion of these transients is out of the 
scope of this paper. There are two fainter X-ray sources, 
CXOU J09555 1.2+694044 (J095551.2 hereafter) and CXOU 
J095550.6+694044 (J095550.6 hereafter) located at about 2 
arcsec to the south of J095551.0 forming a triangle (Fig. 1). 
In contrast to J095551.0, these two sources are always active. 

3.2. X-ray Spectroscopy 

We performed spectral analysis for all Chandra ACIS data 
by using XSPEC vll.3. We also used CIAO's Sherpa for 
independent check. In three of the observations, the X-ray 
sources are located near the aim point and we can use a circu- 
lar extraction region with radii of 0.8- 1.3 arcsec depending 
on the contamination of nearby sources. The relatively small 
extraction radii were used because the three X-ray sources that 
we are interested in are close to each other. For the remaining 
four observations, our targets were off-axis and therefore we 
used an elliptical region to extract the spectra. For the back- 
ground, we used a nearby source free region. We rebinned 
the 0.3-7 keV spectra with at least 20 counts per spectral bin, 
and used \ 2 statistics to find the best-fitting parameters. Cor- 
responding response files were generated using CIAO. 

For the transient (J095551.0), the X-ray spectra can be ad- 
equately fitted with an absorbed power-law model. The spec- 
tral parameters are listed in Table 2. In general, the photon in- 
dex varies between 1.3 and 1.7 while the Nh is about 3 x 10 22 
cm" 2 , consistent with the extinction measured with near IR 
observations (Alonso-Herrero et al. 2003). It is worth noting 
that during observation 9, the source suffers mild pile-up. In 
this case, we included a pile-up model in spectral fit yield- 
ing a pile-up fraction of <~ 15%. For consistency check, we 
also applied a pile-up model for other observations. Pile-up 
affects some of the observations (in particular observations 4, 
6, 7 and 8) with a maximum pile-up fraction of ~ 10%. Our 
spectral fit of observation 10 is different comparing to Kaaret 
et al. (2006). This is likely due to the diffuse background 
and contamination of nearby sources. We used a nearby dif- 
fuse emission region as the background and the fit was per- 
formed with the background subtracted spectrum. We also 
used a smaller extraction region to reduce the contamination 
of nearby sources. In fact, Kaaret et al. (2006) required an ad- 
ditional very soft component to fit the spectrum indicating the 
contribution of diffuse emission and nearby sources. We also 
performed a fit with a bigger extraction region and without 
background subtraction. The result is consistent with Kaaret 
et al. (2006). Since five of the observations (4, 6, 7, 8, and 10) 
are off-axis, the spectra of the transient are contaminated by 
J095551.2 and J095550.6. The contamination is particularly 
serious in observations 4, 6, 7, and 8, and may result the rel- 
atively hard spectra of these observations. To verify the con- 
tamination, we extracted combined spectra of all three sources 
using observations 11 and 12 for which the sources are well 
resolved. While the X-ray flux is dominated by the transient, 
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FIG. 1. — Chandra 0.3-7 keV images of the central 45" X 45" region of M82 as seen on 1999 September 20 (Left; observation 1) and 2002 June 18 (Right; 
observation 9). Both figures have the same spatial scale. The locations of the three ULXs discussed in this paper are marked. We also indicate the position of 
M82 X-l. The images has been slightly smoothed with a 0.5"ct Gaussian function. 



TABLE 2 

Power-law Spectral Model for the ULXs 



Observation 



J095551.0 



J095551.2 



J095550.6 







r 


Lx 


X Z ldof 


Nh" 


r 


Lx 


X'ldof 




r 




x'/dof 


1 

9 
10 
11 
12 


3 63+ - 75 

n in+0.23 
J - ly -0.17 
3 56 +<u9 


1 74+0.51 
1 -'*-0.51 
1 47+0.12 

i <;9+0.09 
l J -0.08 
, 97+0.18 

1 - z '-o.is 


7.8 ±0.2 
12±0.1 
13±0.1 
11 ±0.2 


0.8/69 
1.0/293 
1.1/331 
0.9/126 


,1+5.9 
zl -2.3 

1 -3.2 

1 3+ 2 ' 6 
1J -5.1 
,9+10.5 
lz -46 


2 15+ - 95 

1 ()7+166 
1A "-0.61 

66+ ' 43 
u.oo_ 048 

61 +18 ° 


17 ±0.7 

3.2 ±0.2 

3.6 ±0.09 

3.3 ±0.2 


1.0/81 
0.9/17 

1.4/70 
1.2/16 


2 86+ 045 

9 41 +0.69 
z - H1 -0.60 

2 56+ 035 
Jyy -i 10 


OJ -0.33 
9 81+0.58 
2 - S1 -0.52 

9 79+O.3I 
z -' z -0.28 

4 46+ 112 
^"-0.93 


2.4±0.1 
2.0±0.1 

I. 5 ±0.04 

II. 6±0.6 


0.9/48 
1.1/20 

1.5/58 
1.7/12 



NOTE. — All quoted uncertainties are 90% except for the luminosities which are \u. 

The spectrum of J09555 1 .0 in observation 9 suffered pile-up and a pile-up model was applied during spectral fit. 
" in units of 10 22 cirT 2 

* 0.5-10 keV unabsorbed luminosity in units of 10 39 ergs s" 1 (assuming d=3.6 Mpc). 



the X-ray spectra become significantly harder with a photon 
index of ~ 1 . This indicates that the hard spectra of the three 
contaminated observations are likely due to nearby sources. 
Furthermore, mild pile-up may also affect the spectra. We 
therefore do not include the spectral fits in Table 2. The best 
fitting power-law spectrum of observation 1 1 is shown in Fig- 
ure 2. 

J095551.2 is the next brightest source near the transient. 
We first fitted the spectra with an absorbed power-law model 
and the spectral parameters are shown in Table 2. Three of 
the fits are acceptable and the spectra turn over at about 4 keV 
suggesting very high absorption. The best fitted Nh is about 
( 1 - 2) x 1 23 cm -2 which is an order of magnitude greater than 
the other two nearby sources. In addition, all spectra are very 
hard with T < 1 except for observation 1. For observation 
11, a soft excess is clearly seen in the spectrum (Fig. 2) and 
the fit is much poorer than the others. Indeed, soft excess is 
seen in all spectra but with larger uncertainties due to shorter 
exposure time or smaller collecting area of ACIS-I below 2 
keV. Soft excess is a common feature of AGN with an ion- 
ized absorber. We then refitted the spectrum (observation 11) 
with an additional ionized absorber (absori model in XSPEC; 



Zdziarski et al. 1995). The fit is acceptable with a reduced 
X 2 of 1 .06, and the best-fit photon index steepens to 2 with a 

cm" 2 (see Fig. 3); 



-,23 



ergs cm 2 s 1 . 



large absorbing column of Nh = 9 x 10 
the absorbed 0.5-10 keV flux is 10" 12 i 

Apart from observation 1 1, the X-ray spectra of J095550.6 
can be fitted with an absorbed power-law with N H ~ 3 x 10 22 
cm" 2 and a photon index of > 2.7. Observation 1 1 has the 
longest exposure time and the X-ray spectrum of J095550.6 
is clearly more complicated. It cannot be fitted with any sin- 
gle component model. Instead, a combination of Raymond- 
Smith model and power-law model is required. In addition, 
emission lines are clearly seen in the X-ray spectrum (Fig. 
4). The X-ray spectrum indicates that J095550.6 is either a 
nearby foreground star or a supernova remnant in M82. We 
will show in §4.3 that J095550.6 is indeed a young supernova 
remnant in M82. 

3.3. X-ray Variability 

With the spectral fits, we can estimate the X-ray fluxes of 
the three sources and study the long-term variability. We limit 
our analysis to Chandra data because the three sources as 
well as M82 X-l are not resolved with Einstein, ROSAT, and 
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FIG. 2.— Power-law spectral fits of J095551.0 (top), J095551.2 (middle), 
and J095550.6 (bottom). Spectra are from observation 11. See Table 2 for 
spectral parameters. 

XMM-Newton. The luminous X-ray source, J095551.0, dis- 
plays strong variability on the timescales of months (see Fig. 
5). In particular, the source was not detected in 1999 Septem- 
ber and reappeared in 1999 October. It was below the detec- 
tion limit again in 2000 January and then turned back on in 
2000 March. Figure 5 shows the long-term X-ray lightcurve 
of J095551.0 from 1999 September to 2005 August. For the 
HRC-I observations and observations 4, 6, 7, and 8, we es- 
timated the X-ray flux by assuming an absorbed power-law 
model with Nh = 3 x 10 22 cm" 2 and a photon index of 1.5. 
When the source is active, the X-ray luminosity shows very 
little variability at (7- 13) x 10 39 ergs s" 1 . We determined the 
90% upper limit when the source was undetected; a nearby 
diffuse emission region was used as the background since 
the source is contaminated by strong diffuse emission. We 
note that Feng & Kaaret (2007) reported a much lower upper 
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FIG. 3. — Chandra spectrum (from observation 11) of J095551.2 with an 
ionized absorber plus power-law model (Nh ■ 



■- 9 x 10 23 cnT 2 , T = 2, 



tion parameter £ = 882, \ /dof = 1.06/69). We fixed the absorber tempera- 
ture and Fe abundance at 10 s K and solar values, respectively as the fit was 
not sensitive to these parameters. 
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FIG. 4. — Chandra spectrum (from observation 11) of J095550.6 with 
an absorbed Raymond-Smith plus power-law model (Nh = 2.9 X 10 22 cm -2 , 
kT RS = 0.9 keV, T = 2.01, x 2 /dof = 1.1/56). 

limit (without statistical significance) by assuming the bright- 
est pixel around the source region. 

J09555 1 .2 has a soft spectrum (r = 2) during the first ob- 
servation and it becomes much harder (r < 1) in subsequent 
observations. Because of the high Nh, the flux is very sensi- 
tive to the photon index. For instance, if we fit the spectra with 
an ionized absorber plus power-law model fixing the spectral 
parameters except for the normalization as in observation 11, 
the luminosities are ~ 10 40 ergss -1 . Hence, the source does 
not show significant variability. 

For J095550.6, the apparent softening during the last obser- 
vation is likely an artifact because of the low count rate. We 
used the same Raymond-Smith plus power-law model as in 
observation 1 1 and the spectrum can be fitted equally well. 
The resulting luminosity is about 2.4 x 10 39 ergss -1 , con- 
sistent with other observations. Therefore, the X-ray flux of 
J095550.6 is consistent with being constant. 

We also study the short-term variability of our targets. We 
extracted the source and background lightcurves from the 0.3- 
7 keV event files except that there is no energy filter for HRC- 
I data. We applied similar procedure as discussed in § 3.2 to 
define the source and background regions of our targets. All 
three sources do not show significant variability on timescale 
of hours. We show the short-term lightcurves of the transient 
in Figure 6. 
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FIG . 5 . — Long-term X-ray lightcurve of the ultraluminous X-ray transient, 
J095551.0. The luminosities are determined by spectral fits (see Table 2). 
For the HRC-I observation and observations 4, 6, 7, and 8, we assume an 
absorbed power-law spectral model with N H = 3 X 10 22 cirT 2 and V = 1.5. 
For non-detections, 90% upper limits are shown. 

3.4. Radial Profile 

We investigated the spatial extent of the ULXs using ob- 
servation 11 for which the sources have the highest num- 
ber of counts and are well resolved. Soft (0.3-3 keV) and 
hard (3-7 keV) band counts were extracted from energy fil- 
tered images. We also modelled the point spread functions 
(PSF) of each source using Chandra Ray Tracer (ChaRT) and 
compared them with the measured radial profile of the three 
sources. Counts from these images were extracted in identi- 
cal manner and normalized by the innermost nuclear annulus. 
The resulting radial profiles are shown in Figure 7. 

The surface brightness distributions of all sources but 
J095551.2 in the hard and soft bands are very similar. They 
are centrally peaked and lie above background level out to ra- 
dius ~ 2" with the exception of J09555 1 .2 for which the con- 
tamination from background is significant in the soft band. 
The radial profiles were compared with the PSF models; we 
do not find significant difference based on a Kolmogorov- 
Smirnov (K-S) test for all sources but J095551.2. 

3.5. Near IR Imaging 

In order to compare the Chandra and HST NICMOS im- 
ages, we first aligned the two images. For the X-ray image, we 
used the HRC-I observation taken on 1999 October 28 (obser- 
vation 3) as the reference frame because it has a wide field- 
of-view and moderate exposure. Furthermore, all three tar- 
gets were active during this observation. We used CIAO tool 
wavdetect to detect X-ray sources in the HRC-I image. We 
compared the X-ray source list with the 2MASS catalog, and 
looked for coincidence of bright and isolated stellar objects. 
We found one star (2MASS 09551494+6936143) that is < 1" 
from the corresponding HRC-I position. From the ACIS-I ob- 
servation (observation 9), the X-ray colors of the X-ray emit- 
ting star indicate that it has a very soft X-ray spectrum (84% 
of the source counts come from < 1 keV with no counts above 
2 keV), consistent with a very soft X-ray source (Di Stefano 
& Kong 2004). The X-ray radiation is therefore likely due to 



the coronal emission from a foreground star. The star has a 
R magnitude of 10.1 (Monet et al. 2003). We calculated the 
X-ray to optical flux ratio as \og(fx / fid = log fx + 5 .67 + 0AR 
(Hornschemeier et al. 2001). With a count rate of 9.4 x 10" 4 
c/s in the ACIS-I detector and assuming a Raymond-Smith 
model with kT RS = 0.3 keV and N H = 4x 10 20 cm -2 (the Galac- 
tic value toward the direction of M82), the 0.3-10 keV flux is 
10" 14 ergs cm -2 s" 1 and the corresponding fx/fR is 5.1 x 10" 5 , 
consistent with a foreground star. 

Based on the 2MASS counterpart, the boresight correc- 
tion that needs to be applied to the X-ray source positions 
is 0.87 ±0.56 arcsec in R.A. and 0.73 ±0.27 arcsec in decl.; 
the uncertainties are a quadratic sum of the errors on the X- 
ray and 2MASS positions. To study the IR environment of 
the X-ray sources, we plot on the NICMOS image (Figure 
8) the corrected X-ray positions with error circles given by 
the quadratic sum of the positional uncertainty for the X-ray 
source (0.032" for the transient and 0.063" for the other two 
sources), the uncertainty in the optical astrometry (2MASS to 
ACS astrometry and ACS to NICMOS astrometry; 0.2"), and 
the uncertainty in the X-ray boresight correction (0.62"). In 
the figure, we also plot the locations of several known super- 
star clusters and M82 X-l . As expected, this area shows many 
star forming regions with the presence of super-star clusters. 
The three luminous X-ray sources as well as M82 X-l are 
located near star clusters. In particular, a near IR source is 
at the center of the error circle of J95551.0. The source is 
marginally resolved in the F160W image with a half-light ra- 
dius of 0.6 pc by fitting with a King model (N. McCrady, pri- 
vate communication). At that size the cluster likely has a mass 
of < 10 5 M Q . J095551.2, however, does not seem to have any 
counterpart. While there is no obvious counterpart within the 
X-ray error circle of J095550.6, some unresolved IR emission 
is seen. Furthermore, a super-star cluster, MGG-8, is just 
outside the lcr X-ray error circle locating 1 arcsec from the 
source. It is also worth noting that M82 X-l is located just 
outside the super-star cluster, MGG-1 1. 

4. DISCUSSION 

4.1. J095551.0: An Ultraluminous X-ray Transient 

The most intriguing behavior of J095551.0 is the X-ray 
variability. The source varies from below the detection limit, 
- 2.5 x 10 38 ergs s" 1 in the 0.5-10 keV band, to - 10 40 ergs 
s" 1 on timescales between observations of <~ 2 months. Fur- 
thermore, the source shows recurrent outbursts. Recurrent ul- 
traluminous transients are not common in nearby galaxies. By 
comparing ROSAT and XMM-Newton observations, Winter et 
al. (2006) found that most of the ULXs are persistent sources 
with less than a factor of 3 in flux variation over the timescale 
from ROSAT to XMM-Newton. Nevertheless, ultraluminous 
X-ray transients have been found in NGC 3628 (Strickland 
et al. 2001), M74 (Soria & Kong 2002), NGC 300 (Kong 
& Di Stefano 2003), NGC 253 (Bauer & Pietsch 2005), and 
M101 (Kong et al. 2004; Kong & Di Stefano 2005). Two 
of these sources (NGC 300 and M101) are ultraluminous su- 
persoft sources with kT < 0.1 keV. The ones in NGC 3628 
and M74 are typical power-law sources with T <~ 2 while the 
ULX in NGC 253 can be described with a bremsstrahlung 
model with kT = 2.2 keV. J095551.0, however, has a much 
harder spectrum. These sources also show diverse luminosity 
range. The sources in NGC 253, M74, and NGC 300 only 
reach L x > 10 39 ergs s" 1 suggesting that an IMBH is not nec- 
essary. On the other hand, the sources in NGC 3628 and M 1 1 
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FIG. 6. — Chandra short-term lightcurves of the transient, J095551.0, when it is active. The time resolution of each plot is 500 s. We do not show observation 
8 due to its short exposure time (1.2 ks). We note that the apparent difference in the count rate is due to different detectors and off-axis angle of the source. 



can have 0.3-8 keV luminosities of ~ 5 x 10 40 ergs s -1 with 
bolometric luminosities approaching 10 41 ergs s -1 . While we 
cannot totally rule out a stellar-mass black hole model, a black 
hole of intermediate mass is certainly an attractive scenario. 

A factor of > 50 in flux variation indicates that J09555 1 .0 
is a compact source while recurrent outbursts and hard X- 
ray spectra can rule out the possibility that the source is a 
young supernova remnant. It is also unlikely to be a back- 
ground AGN since AGN normally varies by a factor of < 10 
on timescales of days to months. However, transient AGNs 
are not unusual (e.g., Komossa et al. 2004; Grupe et al. 2004) 
and they belong to a class of AGNs called narrow-line Seyfert 
1 galaxies. These galaxies have X-ray spectra much softer 
(r > 2.5; e.g., Boiler et al. 1996; Grupe et al. 2004) than 
J095551.0. The more likely scenario is that J095551.0 is a 
binary system with a black hole accretor. The high X-ray lu- 
minosity (Lx ~ 10 40 ergs s" 1 ) when it is active indicates that 
it may be an ULX with an IMBH. Assuming the emission is 
isotropic, the X-ray luminosity implies that the compact ob- 



ject is a ~ 100M Q black hole. Many ULXs have a thermal 
component with a temperature of w 0.1 keV which is inter- 
preted as evidence of IMBHs with masses of ~ 100- 1000M Q 
(Miller et al. 2004). J095551.0, however, does not have any 
soft excess and the X-ray spectra are well fitted with an ab- 
sorbed power-law model with photon index T = 1.3 — 1.7. This 
resembles to the low/hard state of Galactic black hole X-ray 
binaries (McClintock & Remillard 2006). It is therefore pos- 
sible that the source can be explained in the framework of the 
advection-dominated accretion flow (ADAF) model. 

More recently, Yuan et al. (2007) apply an ADAF model to 
describe the X-ray emission of M82 X-l and argue that the 
accreting compact object is an IMBH. During the low/hard 
state, the ADAF model predicts that the X-ray luminosity is 
about < 1% - 10% of the Eddington luminosity. If J095551.0 
is accreting at a rate similar to the hard-state of Galactic black 
hole X-ray binaries, this would imply a black hole mass of 
~ 1000- 10 4 M Q . The estimate should be treated with cau- 
tion because we assume that the X-ray spectra of the transient 
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FIG. 7.— The soft (0.3-3 keV; solid points) and hard (3-7 keV; triangles) 
band radial profile of J095551.0 (top), J095551.2 (middle), and J095550.6 
(bottom), compared to a Chandra PSF model (solid curve), from observation 
11. 

ULX are similar to that of Galactic black holes in the low/hard 
state. Indeed, pure power-law spectral model for ULXs is not 
uncommon; Chandra and XMM-Newton observations have 
revealed hard-state ULXs in several nearby galaxies (see e.g., 
Roberts et al. 2004; Winter et al. 2006). Hard-state ULXs 
may be good candidates to IMBHs. If the accreting object is 
instead a stellar-mass black hole in the hard state, the X-ray 



emission must be anisotropic in order to produce such a high 
X-ray luminosity. However, the inner accretion disk of hard- 
state Galactic black hole X-ray binaries are truncated at large 
distances from the black hole and this may be a problem for 
the thick-disk plus central funnels anisotropic radiation model 
(King et al. 2001). Furthermore, the lack of short time vari- 
ability of J095551.0 argues against the relativistic beaming 
model since this would require a very stable jet. It is worth 
noting that the observed photon index of J09555 1 .0 is some- 
time harder than the typical hard-state value of 1 .5 < T < 2. 1 
for Galactic X-ray binaries (McClintock & Remillard 2006). 
It is therefore not clear if we can directly compare with Galac- 
tic black hole X-ray binaries in the hard state. Alternatively, 
it may be a unique state that the ULX is a stellar-mass black 
hole accreting at very high rate. 

In addition to the long-term timing variability, the X-ray 
spectra also vary. Excluding those observations taken with 
off-axis pointings, the photon index is consistent with 1.5- 
1.7 except for the last observation (see Table 2). In the last 
observation, the photon index becomes much harder with 
T = 1.27 ±0.18. Moreover, the spectral change is quite dra- 
matic. The observation taken one day earlier has a photon 
index of 1.52 ±0.10 while the X-ray luminosity does not 
change significantly. On the other hand, the nearby source, 
J09555 1 .2, does not show this dramatic change in the spectra. 
We also check the spectra of a few bright sources in the field 
for the last two observations; none of the sources displays this 
kind of spectral variability. We therefore can conclude that the 
spectral hardening is real. Such an X-ray spectrum is unusual 
for Galactic black hole X-ray binaries. The only exception 
is the fast X-ray nova, SAX J1819. 3-2525 whose spectrum 
is extraordinary hard with T = 0.6-0.9 during a flaring state 
(Markwardt et al. 1999; McClintock & Remillard 2006). The 
last observation of J095551.0is similar to SAX J1819. 3-2525 
but we note that the luminosity is indeed slightly lower than 
that measured one day earlier and the short-term light curve 
does not display strong variability. 

Near IR and radio observations may provide additional 
clues about the nature of J095551.0. From the NICMOS im- 
age (Fig 8), a star cluster is at the center of the X-ray error 
circle of J095551.0, suggesting that the source is associated 
with the cluster. Indeed, young star clusters are ideal places to 
produce IMBHs via the collapse of very massive stars through 
runaway stellar collisions (Portegies Zwart et al. 2004). The 
coincidence of J09555 1 .0 and a star cluster strongly suggests 
that the ULX is produced in the cluster and is consistent with 
a black hole of intermediate mass. Radio emission is also 
detected within the X-ray error circle (Kording et al. 2005; 
Kaaret et al. 2006). The radio source, known as 42.21+59.0, 
has been detected several times (e.g., McDonald et al. 2002). 
The flat spectrum (between 5 and 15 GHz) and extended size 
(4.9 pc) suggest that it is a giant H II region with 94 05 stars 
(McDonald et al. 2002). 

4.2. J095551.2: A Highly Absorbed X-ray Source 

The X-ray spectrum of J095551.2 is very different com- 
pared to that of other nearby sources. The source has an un- 
usually high Nh (> 10 23 cm" 2 ) while the two nearby sources 
as well as M82 X-l (Kaaret et al. 2006) have a N H of 
~3 x 10 22 cm" 1 , consistent with the extinction measured by 
IR observations. Furthermore, J095551.2 has a very flat spec- 
trum (r < 1) with soft excess below 2 keV. The high absorp- 
tion column density may indicate that the source is a back- 
ground AGN. AGN normally has a power-law spectral model 
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FIG. 8.— /OT7NICMOS F160W image of the region around the X-ray sources analyzed in this paper (J095551.0, J095551.2, and J095550.6). The la Chandra 
error circles (0.66 arcsec) are shown. We also label the locations of known super-star clusters and M82 X-l. 



with r ~ 1.7-2 (e.g., Page et al. 2006) while the spectrum 
of J09555 1 .2 is much harder. A hard spectral index of AGN 
may be due to the presence of a reflection component and/or 
a complex absorber (Cappi et al. 2006). In §3.1, we refit- 
ted the spectrum (observation 11) with an ionized absorber 
plus power-law model yielding an acceptable fit. The pho- 
ton index steepens to 2 and the 0.5-10 keV unabsorbed flux 
is significantly higher. The other two spectra (observations 
9 and 12) can also be well fitted with an ionized absorbed 
plus power-law model with the spectral parameters fixed at the 
values determined in observation 1 1 except for the normaliza- 
tion. However, due to the low count rate of these two observa- 
tions, it is not clear if the difference is real. On the other hand, 
the spectrum of the first observation might be different. The 
spectral parameters of observation 1 1 cannot fit the spectrum 
of the first observation. If we only fit all the spectra in the 
range of 2.5-7 keV where soft excess is not crucial, the spec- 
trum of the first observation is much softer (T = 2.9) than the 
others (V < 1). Hence, the spectral change during the first ob- 
servation might be real. Spectral change is also seen in some 
narrow-line Seyfert 1 galaxies such as NGC 405 1 (Guainazzi 
et al. 1996; Ponti et al. 2006; see Leighly 1999 for a re- 
view). Therefore, J095551.2 is likely to be a highly absorbed 
background AGN. It is worth nothing that there is a H II re- 
gion within the X-ray error circle (42.56+580; McDonald et 
al. 2002). 

4.3. J095550.6: A Young X-ray Supernova Remnant 

J095550.6 is unique because the X-ray spectrum cannot be 
fitted with simple spectral models. Instead, a combination of 
a Raymond-Smith model and a power-law model is required 
to describe the spectrum (see Fig. 4). Given the high Nh 



that is consistent with the extinction of nearby region of M82, 
J095550.6 is unlikely to be a foreground star. The X-ray spec- 
trum is not typical for X-ray binaries. We can also rule out a 
background AGN due to its unusual spectrum. The remaining 
possibility is that it is an X-ray luminous supernova remnant. 
In particular, it is evident that the spectrum (Fig. 4) is dom- 
inated by broad emission of Mg XII lines at 1 .4 keV, Si K 
shell lines at 1.8 keV, and S K shell lines at 2.5 keV. Strong 
emission lines are also seen in some luminous X-ray emitting 
supernova such as SN 1978K (Schlegel et al. 2004). Indeed, 
SN 1978K is also the first known supernovae with X-ray lu- 
minosity above 10 39 ergs s" 1 . For direct comparison with SN 
1978K, we refit the spectrum of J095550.6 with an absorbed 
MEKAL + power-law model. The spectral parameters are not 
sensitive to the model, with a best-fit temperature of 0.9 keV. 
This is slightly hotter than SN 1978K (~ 0.7 keV; Schlegel 
et al. 2004). The 0.5-10 keV luminosity of J095550.6 is 
2.2 x 10 39 ergs s _1 . 

Within the X-ray error circle, there is a strong radio source 
known as 41.95+575 (McDonald et al. 2001,2002; Ktirding 
et al. 2005). 41.95+575 is the brightest and most compact ra- 
dio source in M82 and has been detected since 1965 (Muxlow 
et al. 2005). High resolution (3 mas) VLBI imaging shows 
that the source has a double-lobed structure (McDonald et 
al. 2001). The separation of the two brightest components 
is 22.4 mas in 2001 (Muxlow et al. 2005). Furthermore, 
multiple-epoch VLBI observations show that the separation 
is increasing at a rate of 0.24 mas yr" 1 . The radio flux also 
varied over the last 30 years. It has decreased in flux density 
at a rate of ~ 8.8% per year. An age of around 100 years is 
estimated (Muxlow et al. 2005). The simplest explanation 
of the nature of J095550.6 is that it is a supernova event tak- 
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ing place within a high density molecular cloud (McDonald 
et al. 2001). We can also estimate the age of the remnant 
using the X-ray spectral fit. Following Kong et al. (2002), 
assuming J095550.6 is in the adiabatic expansion phase, the 
shock temperature can be written as T s = (0.18 keV)(R / tj,) 2 , 
where R and £3 are the radius (in units of parsecs) and age (in 
units of 1000 yr), respectively. Adopting a radius of 11 mas 
(= 0.19 pc), and T s = 0.9 keV, we obtain t ~ 100 yr which 
is consistent with radio observations. It is worth noting that 
long-term radio observations have suggested that 41.95+575 
may be a radio afterglow of a 100 year old gamma-ray burst 
event (Muxlow et al. 2005). 

5. CONCLUSIONS 

We have used archival Chandra and //ST/NICMOS data to 
study the physical nature of three ULXs near the center of 
M82. We found a recurrent transient ULX, J095551.0 from 
the Chandra data. During its active state, the X-ray lumi- 
nosity is about 7 x 10 39 —1.3 x 10 40 ergs s" 1 and it was turned 
off twice in 1999 and 2000 indicating a factor of > 50 vari- 
ability. This also rules out the possibility that it is a super- 
nova remnant. The X-ray spectra can be fitted with a power- 
law model with photon index T = 1.3-1.7 which is similar 
to Galactic black hole X-ray binaries in the low/hard state. 
We suggest that the X-ray emission might be explained in the 
framework of the ADAF model implying a black hole mass 
of ~ 10 3 - 10 4 M Q . However, we cannot totally rule out that 
the source is in a unique spectral/luminosity state. In partic- 
ular, spectral hardening was seen in one of the observations. 
We also examined near IR images taken with //S77NICMOS. 
We found a star cluster at the center of the X-ray error circle 
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suggesting that the source is associated with the cluster. 

With an unusually high column density Ng > 10 23 cm" 2 and 
a rather flat X-ray spectrum (T < 1) with an ionized absorber, 
it is suggested that the source J09555 1 .2 is likely to be a back- 
ground AGN. The source also shows spectral change similar 
to some narrow-line Seyfert 1 galaxies. In addition, there is 
a H II region known as 42.56+580 associated with the X-ray 
source. 

The source J095550.6 shows an unusual spectrum that can- 
not be fitted with a simple power-law model. Instead, the 
spectrum can be fitted with a Raymond-Smith model, accom- 
panying with broad emissions of Mg, Si, and S, indicating 
that it is a supernovae remnant with an age of ~ 100 years. 
Furthermore, a radio source known as 41.95+575 is found as- 
sociated with J095550.6 within the X-ray error circle. The 
long-term radio observations reveal that 41.95+575 may be a 
radio afterglow of a 100-year old gamma-ray burst event. 
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